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We  discuss  a  dynamic  procedure  that  makes  fractional  derivatives  emerge  in  the  time  asymptotic  limit  of 
non-Poisson  processes.  We  find  that  two-state  fluctuations,  with  an  inverse  power-law  distribution  of  waiting 
times,  finite  first  moment,  and  divergent  second  moment,  namely,  with  the  power  index  /r  in  the  interval 
2</j,<3,  yield  a  generalized  master  equation  equivalent  to  the  sum  of  an  ordinary  Markov  contribution  and 
a  fractional  derivative  term.  We  show  that  the  order  of  the  fractional  derivative  depends  on  the  age  of  the 
process  under  study.  If  the  system  is  infinitely  old,  the  order  of  the  fractional  derivative,  is  given  by  c  =  3 
-fi.  A  brand  new  system  is  characterized  by  the  degree  c  =  /a>-2.  If  the  system  is  prepared  at  time  -ta< 0  and 
the  observation  begins  at  time  t= 0,  we  derive  the  following  scenario.  For  times  0 <t<ta  the  system  is 
satisfactorily  described  by  the  fractional  derivative  with  r  =  3  —  /jl.  Upon  time  increase  the  system  undergoes  a 
rejuvenation  process  that  in  the  time  limit  t^>ta  yields  2.  The  intermediate  time  regime  is  probably 

incompatible  with  a  picture  based  on  fractional  derivatives,  or,  at  least,  with  a  mono-order  fractional  derivative. 

DOI:  10.11 03/PhysRevE. 70.036 1 05  PACS  number(s):  02.50.Ey,  05.40.Fb,  05.60.Cd 


I.  INTRODUCTION 

The  fractional  calculus  has  recently  received  a  great  deal 
of  attention  in  the  physics  literature,  through  the  publication 
of  books  [1,2]  and  review  articles  [3,4],  as  well  as  an  ever 
increasing  number  of  research  papers,  some  of  which  are 
quoted  here  [5-12].  The  blossoming  interest  in  the  fractional 
calculus  is  due,  in  part,  to  the  fact  that  it  provides  one  of  the 
dynamical  foundations  for  fractal  stochastic  processes  [2,4]. 
The  adoption  of  the  fractional  calculus  by  the  physics  com¬ 
munity  was  inhibited  historically  because  there  was  no  clear 
experimental  evidence  for  its  need.  The  disciplines  of  statis¬ 
tical  physics  and  thermodynamics  were  thought  to  be  suffi¬ 
cient  for  describing  complex  physical  phenomena  solely  with 
the  use  and  modifications  of  analytic  functions.  This  view 
was  supported  by  the  successes  of  such  physicists  as  On- 
sager,  who  through  the  use  of  simple  physical  arguments  was 
able  to  relate  apparently  independent  transport  processes  to 
one  another,  even  though  these  processes  are  associated  with 
quite  different  physical  phenomena  [14].  His  general  argu¬ 
ments  rested  on  three  assumptions:  (1)  microscopic  dynam¬ 
ics  have  time-reversal  symmetry;  (2)  fluctuations  of  the  heat 
bath  decay  at  the  same  rate  as  do  macroscopic  deviations 
from  equilibrium;  and  (3)  physical  systems  are  aged.  We 
refer  to  assumption  2  as  the  Onsager  principle  and  show  that 
it  is  tied  up  with  assumption  3. 

Onsager’s  arguments  focused  on  a  system  that  is  in  con¬ 
tact  with  a  heat  bath  sufficiently  long  that  the  bath  has  come 
to  thermal  equilibrium  and  consequently  the  system  is  aged. 
In  statistical  physics  we  know  that  the  bath  is  responsible  for 
both  fluctuations  and  dissipation,  and  if  the  fluctuations  are 
white  the  regression  of  perturbations  of  the  bath  to  equilib¬ 
rium  is  instantaneous.  This  means  that  the  energy  absorbed 
from  the  system  of  interest  by  the  bath,  through  macroscopic 
dissipation,  is  distributed  over  the  bath  degrees  of  freedom 
on  a  very  much  shorter  time  scale  than  the  relaxation  time  of 


the  system.  This  property  is  summarized  in  the  well  known 
fluctuation-dissipation  theorem,  which  has  even  been  gener¬ 
alized  to  the  case  where  the  fluctuations  in  the  bath  do  not 
regress  instantaneously  [15]. 

The  dynamics  of  the  physical  variables  to  which  the  On¬ 
sager  principle  apply  are  described  by  two  different  kinds  of 
equations:  (1)  the  Langevin  equation,  a  stochastic  differential 
equation  for  the  dynamical  variable  and  (2)  the  phase  space 
equation  for  the  probability  density.  Two  distinct  methods 
have  been  developed  to  describe  the  phase  space  evolution  of 
the  probability  density;  the  master  equation  introduced  by 
Pauli  and  the  continuous  time  random  walk  (CTRW)  ap¬ 
proach  of  Montroll  and  Weiss  [16].  The  CTRW  formalism 
describes  a  random  walk  in  which  the  walker  pauses  after 
each  jump  for  a  sojourn  specified  by  a  waiting  time  distribu¬ 
tion  function.  It  was  shown  by  Bedeaux  et  al.  [17]  that  the 
Markov  master  equation  is  equivalent  to  a  CTRW  if  the  wait¬ 
ing  time  distribution  is  Poissonian.  However,  when  the  wait¬ 
ing  time  distribution  is  not  exponential,  the  case  we  consider 
here,  the  equivalence  between  the  two  approaches  is  main¬ 
tained  only  by  generalizing  to  the  non-Markov  master  equa¬ 
tion,  the  so-called  generalized  master  equation  (GME)  [18]. 
Recently,  Metzler  [19]  argued  that  the  GME  unifies  the  frac¬ 
tional  calculus  and  the  CTRW. 

Allegrini  et  al.  [13]  have  shown  that  creating  a  master 
equation  compatible  with  the  Onsager  principle  requires  that 
the  system  be  entangled  with  the  bath  in  such  a  way  as  to 
realize  a  condition  of  stable  thermodynamic  equilibrium. 
This  system-bath  entanglement  is  the  result  of  a  rearrange¬ 
ment  process  that  may  take  an  infinitely  long  time  to  com¬ 
plete,  leading  to  the  replacement  of  the  GME  of  Kenkre  et  al. 
[18],  which  corresponds  to  the  nonstationary  condition,  with 
a  new  GME  compatible  with  the  stationary  condition,  and 
consequently  with  the  Onsager  principle.  Herein,  we  extend 
that  discussion  to  include  the  connection  with  the  fractional 
calculus  in  both  infinitely  and  not  infinitely  aged  condition. 
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This  will  allow  us  to  go  beyond  the  discussion  of  Sokolov 
and  Metzler  [12],  whose  fractional  derivative  refers  to  the 
young  state.  We  shall  find  a  fractional  derivative  operator 
with  a  different  index,  corresponding  to  the  infinitely  aged 
condition,  and  we  shall  show  that  the  condition  of  Sokolov 
and  Metzler  is  an  attractor  for  all  the  systems  that  are  par¬ 
tially  aged  and  not  infinitely  aged. 


Beyond  the  Onsager  principle 

We  approach  the  subject  of  fractional  derivatives  from  a 
perspective  similar  to  that  of  Sokolov  and  Metzler  [12]. 
More  specifically,  we  address  the  problem  of  the  connection 
between  the  GME  [18]  and  the  stationary  version  of  the 
CTRW  [16].  The  GME  considered  by  Allegrini  et  al.  [13]  is 
the  two-site  version  of  the  following  equation: 

~7V(l)  =  -  [  (1) 

at  J0 

where  p(f)  is  the  m-dimensional  population  vector  of  m  sites, 
K  is  a  transition  matrix  between  the  sites,  and  <I>(f)  is  the 
memory  kernel.  The  CTRW  prescription  for  this  process 
yields 

00  r 

P(l)  =  2  dt'if/ll(t')'ir(t-t')M"p(0).  (2) 

n= 0  J  0 

Note  that  i/c„(f)  is  the  probability  that  n  jumps  occurred  and 
that  the  last  jump  took  place  at  time  t=t' ,  implying  the  re¬ 
newal  theory  relation 

<PnU)  =  [  'l'n-l(t-t')<l'\{t')dt' ,  (3) 

do 

where  <//,  (t)  is  the  waiting  time  distribution  function  tp{i) 
introduced  into  the  CTRW  and  i//Q(t)  =  <S(f) .  While  M  is  the 
transition  matrix  connecting  the  sites  after  one  jump  has  oc¬ 
curred,  the  probability  that  no  jump  occurs  in  the  time  inter¬ 
val  (0 ,  f)  is 


V(t)  = 


ip{t')dt' . 


(4) 


The  waiting  time  distribution  function  and  the  memory  ker¬ 
nel  can  be  related  to  one  another  by  taking  the  Laplace  trans¬ 
form  of  the  GME  (1)  and  the  CTRW  (2).  This  comparison, 
after  some  algebra  [13],  yields 

^  uiMu) 

<P(m)K=  (M-l),  (5) 

1  -  i/Au) 

where  I  is  the  m  X  m  unit  matrix  and  the  Laplace  transform 
of  the  function  /(f)  is  /(«).  Here,  as  in  Allegrini  et  al.  [13], 
we  limit  our  discussion  to  the  two-state  case  where 


M  = 


(6) 


and 


K  = 


thereby  reducing  Eq.  (5)  to 


<!>(«) 


utfAu) 
l  -  ijj{a) 


(7) 


(8) 


This  relation  between  the  Laplace  transform  of  the  memory 
kernel  and  the  Laplace  transform  of  the  waiting  time  distri¬ 
bution  function  was  first  obtained  by  Kenkre  et  al.  [18]  and 
is  reviewed  by  Montroll  and  West  [20]. 

In  the  case  when  the  lattice  has  only  two  sites,  a  left  and 
a  right  site,  the  random  walker  corresponds  to  a  dichotomous 
signal  £,  with  the  values  £(?)  =  - 1,  for  the  left  site,  and  £(f) 
=  1,  for  the  right  site.  For  the  sake  of  simplicity,  we  assume 
the  two  states  to  have  the  same  statistical  weight.  Also  in  the 
two-state  CTRW,  if  we  adopt  a  discrete  time  representation, 
the  motion  of  the  random  walker  corresponds  to  a  symbolic 

sequence  {£},  with  the  form  {++  +  +  +  +  H — I - 1  +  +  H - 

- ■••},  which  shows  a  significant  persistence  of  both 

states.  The  waiting  time  distribution  i/At)  is  the  distribution 
of  the  patches  filled  with  either  -i-’s  or  -’s. 

We  assume  symmetry  between  the  two  states  and  a  finite 
first  moment  of  ip(t)  making  it  possible  for  us  to  define  the 
autocorrelation  function  for  the  fluctuations  £(f), 


= 


<mm 

(e) 


(9) 


because  the  process  is  stationary  in  time  [21].  The  problem 
of  how  to  relate  3>^(f)  to  this  persistent  fluctuation  is  delicate 
and  will  be  discussed  with  some  detail  in  Sec.  I  B.  Here  we 
limit  ourselves  to  noticing  that  the  Onsager  principle  by¬ 
passes  the  technical  difficulties  with  the  connection  between 
the  Langevin  and  the  master  equation  pictures,  either  ordi¬ 
nary  or  generalized.  In  fact,  the  master  equation  is  a  prescrip¬ 
tion  to  determine  the  probability  of  occupying  a  given  state. 
In  the  case  under  study  here,  the  master  equation  makes  it 
possible  to  determine  at  time  t,  with  /=  1 ,2.  The  Lange¬ 
vin  equation,  or  in  general  any  theoretical  tool  driving  the 
motion  of  a  variable,  £(t)  in  the  case  under  study  here,  af¬ 
fords  a  convenient  means  to  determine  £(f),  and  conse¬ 
quently  the  autocorrelation  function  of  Eq.  (9). 

In  his  original  work  Onsager  considered  the  case  of  a 
macroscopic  fluctuation  that  regresses  to  equilibrium  through 
the  phenomenological  equations  of  motion.  Here  we  adopt 
an  extension  of  the  Onsager  principle  made  by  Allegrini  et 
al.  [13]  to  the  case  of  two  states,  using  the  probability  of  the 
random  walker  being  in  state  i=  1 ,2,p((f),  at  time  t,  which 
allows  us  to  determine  the  autocorrelation  function  of  Eq.  (9) 
without  any  need  to  establish  the  Langevin-like  picture 
equivalent  to  the  GME  under  study  here.  In  fact,  we  adopt 
the  following  equality: 


<^(f)  = 


Pi(0)-p2(0)’ 


(10) 


which  yields  an  immediate  connection  between  the  GME  of 
Eq.  (1)  and  the  autocorrelation  function  of  Eq.  (9),  provided 
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that  the  system  is  initially  in  the  out  of  equilibrium  state 
corresponding  to  /?  j (0)  — /?2(0)  #0.  Assuming  a  regression  to 
equilibrium  in  such  a  way  as  to  retain  Eq.  (10)  we  obtain 
from  the  GME  (1)  using  the  coupling  matrix  (7) 


The  authors  of  Ref.  [13]  determined  that  the  problem  of 
how  to  make  these  processes  compatible  with  the  Onsager 
principle  could  be  solved  by  expressing  the  CTRW  in  sta¬ 
tionary  form,  resulting  in  the  GME  memory  kernel 


— — -  =  -  2  f  <b(t-t')<S>£t')dt' .  (11) 

dt  J0 

Thus,  the  Laplace  transform  of  the  autocorrelation  function 
can  be  related  to  the  Laplace  transform  of  the  memory  kernel 
by 

$((«)  =  — ■  (12) 

Li  +  2d>(«) 

It  is  interesting  to  notice  that  in  the  Poisson  case,  namely, 
when  ijAj)  is  an  exponential  function  of  time,  the  memory 
kernel  of  the  GME,  given  by  Eq.  (8)  turns  out  to  be  equiva¬ 
lent  to  a  Dirac  S  function  of  time,  thereby  implying  that  the 
bath  responsible  for  the  fluctuations  of  the  variable  £  has  a 
time  scale  infinitely  smaller  than  the  system  of  interest.  In 
this  specific  case,  with  the  help  of  Eq.  (12)  we  see  that  the 
autocorrelation  function  <lL(r)  decays  exponentially  with 
time.  This  is  a  condition  behind  the  ordinary  Onsager  prin¬ 
ciple.  Following  the  authors  of  Ref.  [13]  we  want  to  go  be¬ 
yond  the  ordinary  Onsager  principle. 

The  authors  of  Ref.  [13]  studied  the  case  where  the  auto¬ 
correlation  function  departs  from  the  exponential  relaxation 
and  has  the  the  following  time  asymptotic  property: 

(13) 

The  waiting  time  distribution  ip{t)  corresponding  to  this  au¬ 
tocorrelation  function  has  the  following  time  asymptotic 
property  [13]: 

o~^,  (14) 

with  //>  2  to  fit  the  stationary  condition.  At  first  sight,  one 
might  be  surprised  about  our  decision  to  make  these  complex 
processes  obey  the  Onsager  principle.  Such  processes  have 
exotic  thermodynamical  properties,  and  in  some  cases  they 
are  even  shown  to  be  nonergodic  [22]  and  to  produce  aging 
effects  [23,26].  Another  interesting  effect  emerging  from 
these  processes  was  described  in  Ref.  [24].  These  authors 
used  a  fractional  Fokker-Planck  equation,  which  is  closely 
related  to  the  GME  used  in  this  paper,  and  they  found  that 
the  response  of  their  GME  to  external  perturbation  is  quite 
different  from  the  response  of  the  corresponding  CTRW,  in 
conflict  with  the  fact  that  their  GME  is  equivalent  to  a 
CTRW  in  the  absence  of  perturbation.  All  these  surprising 
properties,  however,  refer  to  the  case  jll<  2,  where  no  invari¬ 
ant  measure  exists.  The  case  /x  >  2 ,  under  study  here,  is  in 
principle  compatible  with  the  Onsager  principle,  and  as  a 
consequence  our  request  is  not  absurd.  Nevertheless,  we 
shall  see  that  the  Onsager  principle  requires  that  the  system 
is  infinitely  aged,  an  ideal  condition,  and  that  an  even  appar¬ 
ently  negligible  departure  from  this  condition  yields  a  strik¬ 
ing  effect:  a  rejuvenation  process. 


«[1  -  (/'(m)] 


$(m)  =  — 

-  2[1  -  ip(n)\  +  m[1  +  < /i(u)\t 

where  t  is  the  average  waiting  time. 


f 

Jo 


tip{t)dt. 


(15) 


(16) 


The  form  of  the  memory  kernel  given  by  Eq.  (15)  is  consis¬ 
tent  with  the  equation  of  motion  for  the  autocorrelation  func¬ 
tion  (11),  and  consequently  Eq.  (15)  is  equivalent  to 


<*>(«)  =  l~( 


1 


2\$f(«) 


(17) 


B.  Theoretical  waiting  time  distribution 

We  have  to  remind  the  reader  that  the  stationary  autocor¬ 
relation  function  of  £  is  not  related  directly  to  i//(  1) .  Zumofen 
and  Klafter  [25]  provided  a  prescription  for  deriving  the  cor¬ 
responding  equilibrium  autocorrelation  function  of  £  from 
i/i(f).  Their  result  rests  on  the  observation  that  ip{t)  is  an 
experimental  function,  evaluated  by  observing  the  time  du¬ 
ration  of  the  two  states.  The  connection  with  renewal  theory 
is  established  by  assuming  that  the  time  duration  of  a  state  is 
determined  by  two  processes;  one  is  the  extraction  of  a  ran¬ 
dom  number  from  a  theoretical  inverse  power-law  distribu¬ 
tion  <//(t),  with  the  same  power  index  /x,  and  the  other  is  a 
coin  tossing  procedure  that  determines  the  sign  of  this  lami¬ 
nar  region.  Thus,  a  given  experimental  sojourn  time  in  one  of 
the  two  states  may  correspond  to  an  arbitrarily  large  number 
of  drawings  and  coin  tossings.  Renewal  theory  is  used  to 
relate  the  autocorrelation  function  ^(f)  to  the  waiting  time 
distribution  function  In  fact,  from  the  renewal  theory 

[21]  we  obtain  the  following  important  result: 

=  7 /  (f'  -t)^*{t')dt' ,  (18) 

where  t  is  the  mean  waiting  time  of  the  -distribution 
density.  It  is  interesting  to  notice  that  this  equation  implies 
that  the  second  derivative  of  the  autocorrelation  function  is 
proportional  to  <//(?), 


£(() 

* 


(19) 


In  Sec.  II  the  departure  point  of  our  calculations  is  given 
by  the  autocorrelation  function  <b ;( 1)  of  Eq.  (18).  In  this  case 
is  convenient  to  assign  to  this  equilibrium  autocorrelation 
function  a  simple  analytical  form.  This  is  done  as  follows. 

First  of  all  we  assign  to  the  following  analytical 

form: 
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=  (20) 

which,  as  can  be  proved  with  the  help  of  Eq.  (24),  is  com¬ 
patible  with  if/l)  getting  the  time  asymptotic  form  of  Eq. 
(14).  This  makes  it  possible  for  us  to  write  t  as  follows: 

r  t* 

T  =  Uf/{t)dt= - -•  (21) 

Jo  M-2 

With  the  choice  of  Eq.  (20)  for  </'  (?),  the  autocorrelation 
function  fT> of  Eq.  (18)  gets  the  attractive  analytical  form 

where 

P=li~  2.  (23) 

Thus,  in  the  case  /x  <  3 ,  the  autocorrelation  function  of  the 
fluctuations  is  not  integrable. 

Zumofen  and  Klafter  [25],  in  addition  to  explaining  with 
clear  physical  arguments  the  connection  between  if/t)  and 
if/*{t),  established  that  the  Laplace  transforms  of  the  two 
functions  are  related  to  one  other  by 


<A*(«) 


2  ifi{u) 

1  +  ij/u) 


(24) 


This  important  relation  allows  us  to  establish  a  connection 
between  r  and  r*,  which  turns  out  to  be 


t=2t.  (25) 

In  conclusion,  there  are  two  different  waiting  time  distri¬ 
butions  with  the  time  asymptotic  behavior  of  Eq.  (14),  the 
experimental  waiting  time  distribution  ifAt)  and  the  theoreti¬ 
cal  waiting  time  distribution  ip*(t).  The  theoretical  distribu¬ 
tion  is  necessary  to  define  the  autocorrelation  function  <lQ(t). 
Thus,  a  theoretical  treatment  involving  the  autocorrelation 
function  will  force  us  to  depend  on  i//(t).  In  this  case  it  is 
convenient  to  adopt  the  analytical  form  of  Eq.  (20)  for 
and  Secs.  II  and  III  will  rest  on  this  choice.  In  other  cases. 
Sec.  IV,  Appendix  A,  and  Appendix  B,  the  theoretical  treat¬ 
ment  will  use  the  experimental  waiting  time  distribution  if/t). 
In  these  cases  it  is  convenient  to  adopt  the  analytical  form 

i/At)  =  ■  (26) 


However,  whatever  choice  is  made,  either  the  analytical  form 
of  Eq.  (20)  for  if/  (t)  or  the  analytical  form  of  Eq.  (26)  for 
if/t),  in  both  cases,  thanks  to  Eq.  (24),  the  two  waiting  time 
distributions  maintain  the  same  time  asymptotic  behavior, 
with  the  same  /x.  So  do  the  two  different  expression  for  the 
equilibrium  autocorrelation  function,  both  decaying  as  Eq. 
(13).  For  the  main  purpose  of  this  paper  the  time  asymptotic 
behavior  is  in  fact  the  property  that  matters. 


C.  Purpose  of  this  paper 

The  purpose  of  this  paper  is  to  prove  that  a  dichotomous 
fluctuation  £,  with  the  waiting  time  distribution  of  Eq.  (14), 
can  be  described  by  a  GME  that,  in  turn,  is  well  represented 
by  a  fractional  derivative  operator.  The  fractional  index  cor¬ 
responds  to  the  realization  of  the  Onsager  condition,  as  sur¬ 
prising  as  this  fact  might  be  in  the  condition  3  >  /x  >  2,  which 
is  so  far  from  the  Poisson  limit  /x=°°.  The  Onsager  principle 
is  not  totally  foreign  to  this  anomalous  physical  condition, 
thanks  to  the  fact  that  /x  >  2  is  compatible  with  the  existence 
of  thermodynamical  equilibrium.  However,  the  time  neces¬ 
sary  to  reach  this  thermodynamical  condition  is  infinite,  and 
the  system,  observed  at  finite  times,  no  matter  how  long, 
shows  a  surprising  rejuvenation  effect.  Through  this  rejuve¬ 
nation  effect,  the  fractional  order  compatible  with  the  On¬ 
sager  principle  slowly  turns  into  that  established  by  Sokolov 
and  Metzler  [12], 

The  ouline  of  the  paper  is  as  follows.  In  Sec.  II,  using  the 
inverse  Laplace  transform  of  Eq.  (17)  we  determine  the  un¬ 
known  memory  kernel  <T(f),  making  it  possible  to  discuss 
how  to  express  the  GME  in  terms  of  fractional  derivatives. 
The  case  where  2</x<3  is  compared  to  the  recent  work  of 
Sokolov  and  Metzler  [12].  We  find  that  the  index  of  the 
fractional  derivative  is  3  -/x,  rather  than  /x- 2,  as  predicted 
by  Sokolov  and  Metzler.  In  Sec.  Ill,  we  prove  that  this  dif¬ 
ference  in  index  is  due  to  the  fact  that  we  adopt  a  stationary 
condition,  while  Sokolov  and  Metzler  do  not.  In  Sec.  IV,  we 
also  prove  that  in  the  case  of  a  finite,  rather  than  infinite,  age 
our  GME  makes  a  transition  from  the  (3-/x)th  to  the  (/x 
-2)th  order.  The  stationary  case  becomes  stable  only  in  the 
limiting  case  of  infinite  age.  In  Appendix  A,  the  interested 
reader  can  find  details  on  how  to  establish  the  order  of  the 
fractional  operator  in  the  GME,  in  the  whole  range  1  </x<  3, 
when  only  the  brand  new  condition  of  Eq.  (8)  is  considered. 
The  accuracy  of  the  analytical  expressions  that  we  use  in 
Sec.  IV  to  illustrate  the  rejuvenation  process  is  discussed  in 
Appendix  B. 

In  conclusion,  on  the  one  hand  we  shed  light  on  the  mean¬ 
ing  of  the  work  of  Allegrini  et  al.  [13],  which  is  proven  to  be 
a  subordination  to  a  Markov  master  equation  through  the 
stationary  distribution  of  first  exit  times.  On  the  other  hand, 
we  extend  the  approach  to  systems  of  any  age  and  reveal  the 
phenomenon  of  a  continuous  time  random  walk  with  rejuve¬ 
nation.  To  accomplish  this  dual  role  we  rely  heavily  on  the 
results  recently  obtained  by  Barkai  [26]  and,  to  a  lesser  ex¬ 
tent  the  results  of  Allegrini  et  al.  [13].  However,  this  allows 
us  to  reveal  a  dependence  of  the  fractional  derivative  order 
on  the  aging  and  rejuvenation  process,  which  was  not  previ¬ 
ously  identified. 


H.  THE  INVERSE  LAPLACE  TRANSFORM  OF  THE 
MEMORY  KERNEL 

To  establish  the  form  of  the  unknown  memory  kernel 
<I>(f),  we  make  a  few  preliminary  observations.  First  of  all, 
we  note  that  through  Eq.  (17)  we  establish  a  direct  connec¬ 
tion  with  the  autocorrelation  function  <lL-(f)  and  that  this 
auto-correlation  function  is,  in  turn,  directly  related  to  the 
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waiting  time  distribution  through  Eq.  (18).  Thus,  with 
no  loss  of  generality  for  the  reasons  illustrated  in  Sec.  I  it  is 
convenient  to  refer  to  i//(t)  rather  than  to  ifii t)  and  conse¬ 
quently,  according  to  the  prescriptions  illustrated  at  the  end 
of  Sec.  I  B,  to  the  analytical  form  of  Eq.  (20). 

For  simplicity,  we  set  T*  =  1  throughout  this  section.  Thus, 
the  Laplace  transform  of  the  autocorrelation  function  is  [27] 

r(^il/V-E“.i),  (27) 

where  0  <  /3  <  1 ,  given  the  fact  that  we  are  considering 
2</i<3,  and  Elp_ j  is  a  generalized  exponential  function  [2]. 
Thus,  <&£(«)  diverges  as  m— > 0  and  Eq.  (12)  yields  <f>(0)  =  0. 
We  explore  the  opposite  limit  u  -  * °°  using  Eq.  (15),  which 
yields  <&(«)=  1  /r=  1/2 t.  In  the  time  representation,  the  lat¬ 
ter  limit  is  equivalent  to  <i>(f)~  S(t)/2r*  for  f— >0.  Therefore, 
we  segment  the  Laplace  transform  of  the  GME  memory  ker¬ 
nel  into  two  parts  as  follows: 

Mu)  =  A:  +  •  (28) 

2  T 

The  first  term  models  the  short-time  limit,  while  the  second 
term  is  responsible  for  the  long-time  behavior.  In  the  time 
representation  we  have 

d>(t)  =  ^  +  <LaU).  (29) 

2  T 

Note  that  this  division  of  the  memory  kernel  into  a  white- 
noise  contribution  and  a  slow  term  corresponds  to  a  similar 
partition  made  by  Fuliriski  [28]. 

Thus,  for  the  time  evolution  equation  of  the  autocorrela¬ 
tion  function  of  £(f),  we  derive  the  following  equation: 

=  _  Icp  (f)  _  2  f  c Pa(t  -  .  (30) 

dt  t  J0 

Using  Eq.  (9)  and  substituting  into  it  the  explicit  expression 
of  t*  as  a  function  of  //,  after  some  algebra,  we  obtain 

/  T  Y  /3t  ^  fit  ['  . 

\t  +  T/  T(t  +  T)~  ^ T(t+T)~~2JQ 


FIG.  1.  The  slow  component  of  the  memory  kernel  <I>(t),<ha(f), 
as  a  function  of  time.  The  black  dots  denote  the  result  of  the  nu¬ 
merical  inversion  of  the  expression  in  the  Laplace  transform  result¬ 
ing  from  Eqs.  (15)  and  (28)  for  /3=0.5;  the  continuous  line  is  the 
analytical  approximation  given  by  Eq.  (36). 


In  the  limiting  case  f— °  we  neglect  the  second  term  on 
the  left-hand  side  of  this  equation.  This  is  a  natural  conse¬ 
quence  of  the  assumption  that  the  memory  kernel  must  tend 
to  zero  with  a  negative  tail  as  an  inverse  power  law.  With  this 
assumption  it  is  straightforward  to  prove  that  the  modulus  of 
the  first  term  becomes  much  larger  than  that  of  the  second 
term  on  the  left-hand  side  of  this  equation.  The  consequences 
of  this  crucial  assumption  are  supported  by  the  numerical 
results  depicted  in  Fig.  1.  With  this  assumption  Eq.  (32)  sim¬ 
plifies  to 

-2/„' 03) 


which  can  be  solved  by  means  of  the  fractional  calculus  [2]. 
We  use  the  Riemann-Liouville  (RL)  definition  of  the  frac¬ 
tional  integral: 


D;q[f(t)] 


1  r 

r^JoU-O1-9’ 


(34) 


which  is  the  antiderivative  of  the  fractional  derivative  with 
order  q,  with  q<  1.  Consequently,  for  /3 <  1  we  can  express 
Eq.  (33)  in  terms  of  the  RL  fractional  integral 


(3D 

The  two  terms  on  the  left-hand  side  Eq.  (31)  are  positive. 
Due  to  the  negative  sign  on  the  right-hand  side  of  this  equa¬ 
tion  we  conclude  that  it  might  well  be  that  <I>„(f)  is  always 
negative. 

Let  us  concentrate  on  the  case  f3<  1:  using  the  autocorre¬ 
lation  function  (f > f( t)  of  Eq.  (22)  (with  T  =  1)  and  using  the 
change  of  time  variable  t'  + 1  — »f',  we  rewrite  Eq.  (31)  in  the 
form 

ft+1  if1  1 

-2  <Da(i+l-t')7/  +  2  <S>a(t+l-t')—dt' 

Jo  '  Jo  ' 


Pt 

(f+  l)P+l' 


(32) 


1  Pt 

2T(1  -  /3)(t+  1)^+1  ’ 


so  that  inverting  this  equation  we  have  the  formal  expression 
for  the  slow  part  of  the  memory  kernel 


®a(*)  =  - 


L 

1 

2T(1  -  fj) 


D 


1-/3 


Pt 

(t+l)P+1 


D)~\pn J»fi/8fl(f)]. 


(35) 


We  denote  the  autocorrelation  function  with  the  inverse 
power  (/3+1)  by  the  symbol  (b^/j+iD-  Carrying  out  the  re¬ 
quired  calculations  on  the  right-hand  side  of  Eq.  (35),  we 
obtain  (see,  for  example.  West  et  al.  [[2],  p.  90], 
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$„(*) : 


P 


2Y(1  -  (3)Y(1  + /3)  {t  +  l)2’ 


P<  1-  06) 


A  comparison  with  numerical  inversion  of  the  kernel  is 
shown  in  Fig.  1. 

For  the  sake  of  completeness,  it  is  worth  noticing  that  we 
can  proceed  in  a  similar  way  in  the  case  /3>  1  also.  In  this 
case  we  find,  for  the  contribution  <I>„(f)  of  the  GME,  the 
following  time  asymptotic  behavior: 


=  - 


PiP~0 


u+  iy*+ 


/?>!. 


(37) 


HI.  THE  EMERGENCE  OF  FRACTIONAL 
OPERATORS 

In  this  section  we  show  that  in  the  two-site  case  we  are 
discussing  the  GME  has  the  form  of  a  transport  equation, 
with  two  terms  on  the  right-hand  side.  The  first  has  the  form 
afforded  by  the  ordinary  master  equation  and  consequently 
satisfies  the  Onsager  principle,  giving  a  relaxation  dependent 
on  the  average  waiting  time  of  the  CTRW.  The  second  term 
corresponds  to  a  fractional  derivative  in  time,  and  extends 
the  Onsager  principle  to  the  case  of  a  relaxation  with  a  fat 
tail.  To  obtain  these  results,  we  use  what  we  have  learned  in 
the  preceding  section. 

First  of  all,  since  we  are  dealing  with  the  two-site  case, 
using  the  form  of  M  and  K  matrices  with  j=x=-l  and  y 
=  1,  we  rewrite  Eq.  (1)  in  the  following  form: 


f  $(l- T)\PiiT) -Piif)]dr, 
0 

(38) 

=  J 

f  $(!- PiPiif) -p2{f))dT. 
0 

(39) 

The  memory  kernel  $  is  related  to  the  autocorrelation  func¬ 
tion  of  the  dichotomous  variable  <I>^  through  Eq.  (17).  Insert¬ 
ing  Eq.  (17)  into  the  Laplace  transform  of  the  set  of  the 
two-site  dynamical  equations,  solving  the  resulting  set  of 
equations,  and  taking  the  corresponding  inverse  Laplace 
transforms  yields  the  solutions 

{1  -<Ff(f)[/72(0)-pi(0)]}, 

(40) 

[1  +  <I)f(f)(p2(0)-pi(0))]. 

(41) 

2  T 


(42) 


with  t*=T*/(/u-2)=1//3,  thanks  to  the  fact  that  we  set  T* 
=  1 .  Substituting  the  decomposition  of  the  memory  kernel 
into  Eq.  (38),  we  obtain 


dpi(t) 


dt 


It 


S' 

Jo 


®a(t- r)lP2ir) -Pi{r)]dT. 


(43) 


Writing  <I/(f)  as  the  derivative  of  an  as  yet  unspecified  func¬ 
tion  /(f),  and  using  the  property 


£> 


—  I  fU  -  T)g(T)dT  =  f(0)g(t)  + 


with  f  (t)  =  id/ dt)f(t),  we  obtain 
dp  i{t)  p2(t)~Pi(t ) 


fr 

Jo 


(t-  T)g(r)dT 


(44) 


dt 


2  T 

X[pi(r)  - pi(r)\dT. 


f(0)[p2(t)  - P\(t)]  +  —  I  fit ~  t) 


-f 

tj  0 


(45) 


We  also  found  that  in  the  case  0  </3<l  the  asymptotic  be¬ 
havior  of  the  memory  kernel  is  expressed  by 

.  M  P  f3 


2T{\- p)Y{\  + 0){t+\)2 


2Y(l-(3)Y(l  +  f3) 


fP  2  t  — >  OO 


1. 


(46) 


1 


tP- 1 


(47) 


Rewriting  Eq.  (47)  as 

$  (t)  =  - - 

2Y{\-  p)Y{p){p-\)df 

we  identify  /(f)  with  [1/2(1  -/3)r(l -/3)r(/3)]f^-1  for  f— >oo. 
Choosing  /(0)  =  0  and  using  the  properties  of  the  F  function, 
we  assign  to  the  time  asymptotic  equation  of  motion  the 
form 


dp  lit)  Piit)-p2{t) 


1 


dt  2 t  2Y(2  -  P)Y{p)  dt 

~  ’dP~l{p2{'t)  ~  Pii^dr 
or,  in  terms  of  the  RL  fractional  integral  (34), 

dp  lit)  Piit)-p2{t)  1 


/' 

Jo 


it 


(48) 


dt 


o  * 

2  r 


2r(2  -  0) 


D)  \p\it)-p2{t)]. 


(49) 


Note  that  these  solutions  can  be  combined  to  yield  the  gen¬ 
eralized  Onsager  principle  given  by  Eq.  (10)  in  terms  of  the 
difference  in  the  probabilities. 

We  now  want  to  find  a  formal  equation  of  evolution  for 
the  probabilities  involving  fractional  operators.  We  know 
from  the  preceding  section  that 


The  same  procedure  applied  to  the  equation  of  motion  for 
Piit)  yields 


dp2jt) 

dt 


P lit) -p lit) 
~  * 

2  r 


1 


2f(2  -  0) 


d)  ^[piit)  -  p2it)\. 


(50) 
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The  difference  between  these  two  equations  yields 


d_ 

dt 


[Pi(t)-p2(t)]  =  - 


P\{t)-p2{t) 

* 

T 


l 

r(2  -  p) 


D^\p  i(t) 


(51) 


clearly  showing  the  two  kinds  of  contribution  to  the  gener¬ 
alized  Onsager  principle.  The  first  term  gives  the  relaxation 
of  the  perturbation  away  from  equilibrium  at  the  macro¬ 
scopic  rate  required  by  Onsager.  The  second  term  gives  the 
additional  slow  relaxation  in  the  form  of  the  fractional  inte¬ 
gral. 


IV.  AGING  ORDER 

We  have  to  remark  that  the  condition  of  Eq.  (12)  refers  to 
the  stationary  condition  explicitly  considered  by  Klafter  and 
Zumofen  [29].  In  this  section  we  prove  that  there  is  a  con¬ 
nection  between  a  system’s  age  and  the  order  of  the  frac¬ 
tional  derivative  in  the  relaxation  process.  A  sign  of  the  de¬ 
pendence  of  the  fractional  derivative  order  on  age  is  given  by 
the  discrepancy  between  the  results  of  Sec.  Ill  and  Ref.  [12]. 
Let  us  compare  Eq.  (49)  to  Eq.  (16)  of  Ref.  [12].  We  obtain 
the  fractional  index  1-/3  rather  than  /3  as  in  the  work  of 
Sokolov  and  Metzler.  In  Appendix  A  we  prove  that  our  time 
asymptotic  approach  to  fractional  derivatives,  in  the  nonsta¬ 
tionary  case  studied  by  Sokolov  and  Metzler,  yields  the  same 
index  as  they  obtain  [12].  Thus,  the  discrepancy  between  our 
prediction  and  the  prediction  of  Sokolov  and  Metzler  de¬ 
pends  on  the  fact  that  we  consider  a  condition  consistent 
with  the  Onsager  principle,  whereas  Sokolov  and  Metzler  do 
not.  Furthermore,  if  the  system  is  not  infinitely  aged,  a  sort 
of  rejuvenation  process  is  expected  to  take  place  that  will 
lead  to  the  fractional  order  of  Sokolov  and  Metzler. 

To  support  our  remarks  concerning  the  relation  between 
aging  and  the  order  of  the  fractional  operator,  here  we  dis¬ 
cuss  how  to  define  a  waiting  time  distribution  of  any  age. 
The  authors  of  Ref.  [13]  have  shown  that  the  waiting  time 
distribution  ip(t)  of  Eq.  (14)  is  obtained  from  the  following 
dynamic  model.  A  particle  moves  in  an  the  interval  I 
=  [0,1]  driven  by  the  equation  of  motion 


at  t  =  -ta+Tx  and  ending  at  t=—ta+  ~]  +  r2,  and  so  on.  The 
waiting  time  distribution  of  age  ta,  denoted  by  iji,  (f),  is  de¬ 
termined  by  the  first  of  these  time  intervals  overlapping  with 
t>  0.  The  time  length  of  that  overlap  is  the  time  length 
whose  distribution  determines  tp,  (t).  We  make  the  assump¬ 
tion  that  the  beginning  of  the  first  time  interval  overlapping 
with  t  >  0  occurs  with  equal  probability  at  any  point  between 
t=-ta  and  /  =  ().  The  validity  of  this  assumption  is  discussed 
in  Appendix  B,  which  establishes  that  this  assumption  is  very 
good  for  ta  — * 0  and  la— In  between  the  asymptotic  limits 
the  resulting  prediction  is  not  exact.  However,  since  it  yields 
simple  analytical  formulas,  we  adopt  this  simplifying  as¬ 
sumption  for  any  age.  Thus,  we  have  that 

J  <Kt  +  y)dy 

'  ,53) 

where  g(ta)  is  the  normalization  factor  defined  by 

g(tj  -  f'ww,  (54) 

do 

and  4'  ( 1)  is  the  probability  that  no  event  occurs  throughout 
the  time  interval  of  length  t.  Using  for  ijii  1) ,  according  to  the 
prescription  adopted  in  this  paper,  the  analytical  form  of  Eq. 
(26),  it  is  easy  to  prove  that  Eq.  (53)  can  be  written  in  the 
form 


(t+T^-it+T+t^ 
^(t)  =  {ji-  2) - ^)_(fa  +  r)(2-rf - 


(55) 


This  formula  proves  that  for  t<ta  the  index  of  the  distribu¬ 
tion  is  fi-1,  whereas  for  ti>ta  the  index  becomes  ji.  This 
result  for  the  age -dependent  waiting  time  distribution  func¬ 
tion  agrees  with  the  predictions  by  Barkai  [26]  and  by  the 
authors  of  Ref.  [13].  Notice  that  the  formula  Eq.  (55)  is 
equivalent  to  drawing  the  initial  condition  for  y  from  an  aged 
distribution  of  this  variable. 

Here,  we  are  in  a  position  to  evaluate  the  waiting  time 
index  at  a  generic  time,  where  we  write  ip,  (f)  as 


dy/dt  =  ayz,  (52) 

with  z>l.  When  the  particle  reaches  the  border  y=l,  it  is 
injected  back  to  an  initial  condition  between  y=0  and  y=l 
with  uniform  probability.  The  age  of  the  CTRW  is  deter¬ 
mined  by  the  distribution  of  first  exit  times.  The  ordinary 
CTRW  is  based  on  identifying  this  distribution  with  ip(t). 
This  means  that  the  CTRW  is  equivalent  to  assuming  that  the 
system  is  prepared  in  a  flat  distribution  at  t= 0,  which  coin¬ 
cides  with  the  beginning  of  the  observation  process. 

Let  us  discuss  now  the  consequence  of  beginning  the  ob¬ 
servation  a  significant  time  after  the  preparation.  Let  us 
imagine  that  the  system  is  prepared  in  a  flat  distribution  at  a 
time  r=-f„<0,  and  that  the  observation  begins  at  t= 0.  This 
means  that  the  flat  distribution  begins  producing  a  sequence 
of  time  intervals  of  size  r,  according  to  the  distribution  of 
Eq.  (14);  more  precisely,  the  time  interval  r,  beginning  at  t 
=  -ta  and  ending  at  t=-ta+T\,  the  time  interval  t2  beginning 


*'St)  = 


MT,ta ) 

(t  +  TjArffW 


(56) 


Using  Eq.  (55)  we  arrive  at  the  following  formula  for  the 
time  dependence  of  the  effective  power-law  index: 


P'e  ffO)  — 


In  [(f  +  r)(1~^  -(t+T+  ra)(1~^] 
ln[r  +  T] 


(57) 


Figure  2  illustrates  the  regression  of  the  effective  power- 
law  index  to  /x  with  two  different  ages,  and  shows  clearly 
that  the  regression  is  slower  for  older  systems.  This  formula 
does  more  than  explain  the  discrepancy  between  Eq.  (49) 
and  Eq.  (16)  of  Ref.  [12].  In  fact,  it  shows  that  it  is  possible 
to  build  a  GME  that  at  short  times  follows  the  prescription 
of  our  GME  and  at  long  times  moves  into  the  basin  of  at¬ 
traction  of  Sokolov  and  Metzler.  This  is  certainly  the  case  if 

t  >*—00 
la 
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FIG.  2.  The  effective  power  index  (/)  as  a  function  of  time, 
for  /r= 2.3.  The  curves  refer,  from  top  to  bottom,  to  fa=100,  1000, 
10  000. 


FIG.  3.  The  Laplace  transform  of  the  ta-old  memory  kernel 
<T,  ( u ),  as  a  function  of  u,  for  /r  =  2.3.  The  curves  refer,  from 

top  to  bottom,  to  ta= 0,0. 1,1,10,°°. 


This  aspect  is  important  and  needs  a  more  exhaustive  il¬ 
lustration.  We  note  that  the  approach  of  Ref.  [13]  can  be 
easily  extended  to  the  case  where  the  distribution  of  first  exit 
times  has  a  finite  age.  It  is  enough  to  follow  the  procedure  of 
Ref.  [13]  and  to  replace  the  first  exit  time  distribution  with 
(A,  (f)  rather  than  with  iAx,(f)  =  <p,  =-dj),  as  done  in  Ref.  [13]. 
The  result  of  this  procedure  yields  for  the  GME  the  follow¬ 
ing  form  for  the  Laplace  transform  of  the  memory  kernel: 

„  mb.  ( u ) 

<M«)  = - .  ‘  . - ■  (58) 

1  +  ip{u)  -  2ipt{u) 

It  is  straightforward  to  prove  that  for  tu= 0  Eq.  (58)  reduces 
to 


&t=  o(w)  = 


Uljj{u) 

1  -  |/i(m) 


(59) 


In  fact,  the  general  prescription  of  Eq.  (53)  immediately 
yields  ip,  (t)  =  i/Xt)  for  ta=0.  To  derive  the  memory  kernel 
corresponding  to  the  infinitely  aged  condition  of  Eq.  (15)  we 
have  to  notice  first  that  Eq.  (53)  yields,  in  accordance  with 
Refs.  [26]  and  [13], 


1  f 

tpJj)  =  ip,  =oo  (f)  =  -  dt'ip(t').  (60) 

"  Tj, 


To  illustrate  the  change  of  the  memory  kernel  with  time, 
notice  that  the  Laplace  transform  of  Eq.  (55)  is 


'V(u){\  -  eL"a)  +  eu>aJ  ‘  e~uy^{y)dy 

fpt  (m)  = - 77V"" - »  (61) 


gUa) 


that  is. 


[1  -  (/<«)](!  -  e~ut“)  +  uelt,“  e~uy^{y)dy 


<A,  (m)  = 


J  0 


ug(ta) 


(62) 


Let  us  substitute  Eq.  (62)  into  Eq.  (58),  from  which  we 
obtain  the  Laplace  transform  of  the  memory  kernel  for  an 
arbitrary  age.  In  Fig.  3  we  show  the  Laplace  transform  of  the 
memory  kernel  corresponding  to  a  number  of  different  ages. 
Moving  from  the  top  to  the  bottom  curve  the  age  increases 
from  the  brand  new  kernel  (ta= 0)  to  the  infinitely  aged  or 
stationary  kernel  (fa=°°). 


V.  CONCLUDING  REMARKS 

The  adoption  of  a  two-state  model  to  generate  anomalous 
diffusion  is  not  unusual  in  the  random  walk  literature.  For 
instance,  Shushin  [30]  also  generated  anomalous  diffusion  by 
means  of  a  two-state  model.  However,  to  properly  locate  in 
the  literature  of  anomalous  diffusion  the  model  of  the  present 
paper,  we  have  to  point  out  that  the  diffusion  generator  in 
Ref.  [30]  is  a  two-state  Markov  equation  modulated  by  time- 
dependent  transition  parameters  that  obey  non-Poisson  statis¬ 
tics.  The  model  studied  in  the  present  paper  might  be 
adopted  to  generate  anomalous  diffusion  as  well,  but  is  quite 
different  from  that  of  Sushin.  The  extended  time  of  sojourn 
in  one  of  the  two  states  would  produce  uniform  motion,  with 
no  randomness  involved,  as  discussed  by  the  authors  of  Ref. 
[12],  which  rests  in  fact  on  the  same  model  as  that  adopted 
here. 

In  this  paper  we  establish  that  an  infinitely  old  system, 
with  the  power-law  index  in  the  interval  2  /jl  =£  3,  yields  the 
fractional  order  1  —  /?=  3  —  /x.  The  prediction  of  Ref.  [12],  on 
the  other  hand,  yields  the  fractional  order  /jl- 2,  correspond¬ 
ing  to  the  brand  new  condition.  If  the  system  is  not  infinitely 
aged,  namely,  the  short-time  behavior  of  the  system 

for  t<ta  is  expected  to  be  that  of  an  aged  system.  At  large 
observation  times,  t>ta,  rejuvenation  begins.  This  can  be 
explained  using  the  dynamical  model  of  Eq.  (52).  In  fact, 
aging  has  to  do  with  the  slow  regression  to  equilibrium,  if  it 
exists,  of  the  variable  y,  which  mimics  a  real  bath  slowly 
regressing  to  equilibrium.  The  aging  effects  discussed  by 
Barkai  [26]  correspond  to  z>2  and  thus  to  fx<2.  In  this 
case  the  dynamic  model  under  discussion  does  not  have  an 
invariant  distribution,  and,  consequently,  any  observation  is 


036105-8 


AGING  AND  REJUVENATION  WITH  FRACTIONAL... 


PHYSICAL  REVIEW  E  70,  036105  (2004) 


done  while  the  bath  is  drifting  toward  a  condition  that  will 
never  be  reached.  This  aging  effect  affects  the  form  of  the 
first  exit  time  distribution,  whose  index  is  /i- 1  rather  than  /jl. 
However,  after  the  first  exit  the  trajectories  are  injected  back 
with  a  uniform  probability,  and  thus  all  the  ensuing  jumps 
are  determined  by  the  ordinary  waiting  time  distribution  < p(t). 

It  would  be  desirable  to  have  an  equation  of  motion  with 
a  fractional  operator  order  that  changes  as  a  function  of  time 
from  3- /jl  to  /jl— 2.  However,  there  are  technical  and  concep¬ 
tual  difficulties  that  make  it  difficult,  if  not  impossible,  to 
realize  this  goal.  In  fact,  according  to  the  perspective  adopted 
in  this  paper,  the  order  of  the  fractional  operator  is  estab¬ 
lished  using  time  asymptotic  arguments.  Thus,  if  ffl<°°  we 
can  associate  a  time  t  <  ta  with  the  order  r= /jl- 2  of  the  frac¬ 
tional  derivative.  This  is  so  because  for  t  <Sfa  the  calculations 
would  be  virtually  equivalent  to  those  done  in  Appendix  A.  It 
is  not  clear  how  to  proceed  when  t  is  of  the  order  of  ta,  this 
being  the  first  reason  why  assigning  a  fractional  derivative 
order  to  any  time  might  be  difficult.  There  also  exists  a 
physical  reason  that  might  make  it  impossible  to  move  from 
the  order  3  -  /jl  to  /jl- 2.  Physically,  this  extremely  extended 
transition  process  might  involve  a  mixture  of  fractional  de¬ 
rivatives  of  different  orders.  It  has  been  shown  [31]  that  the 
Levy  walk  does  not  have  well-defined  scaling,  due  to  aging 
effects.  Similarly,  the  adoption  of  a  fractional  time  deriva¬ 
tive,  with  time-dependent  order,  might  be  inadequate  to  ex¬ 
plore  the  regime  of  transition  from  the  order  3  —  /i  to  the 
order  /i-2.  In  conclusion,  the  fractional  operator  and  its  or¬ 
der  reflect  a  stable  condition,  of  a  brand  new  or  infinitely  old 
system.  The  regime  of  transition  from  the  dynamic  to  any  of 
these  two  thermodynamic  regimes,  and  the  regime  of  transi¬ 
tion  from  the  earlier  to  the  latter  thermodynamic  regime,  is 
not  yet  a  fully  understood  physical  condition,  an  issue  calling 
for  further  investigation. 

It  is  interesting  to  notice  that,  even  if  we  select  /jl>2  and 
consequently  adopt  a  condition  compatible  with  the  station¬ 
ary  condition,  the  effective  index  of  the  first  exit  distribution 
is  located  in  the  nonstationary  region  if  f3<l.  This  is  prob¬ 
ably  the  reason  why  the  memory  kernel  seems  to  share  the 
same  properties  as  those  adopted  in  Refs.  [6,10,32]  to  pro¬ 
duce  subdiffusion.  Notice  that  the  baths  used  by  Lutz  [6]  and 
Pottier  [32]  have  properties  quite  different  from  the  subordi¬ 
nation  perspective  of  Ref.  [10],  even  though  the  relaxation 
process  stemming  from  subordination  [10]  is  quite  similar  to 
that  produced  by  the  non-Ohmic  baths  of  Lutz  and  Pottier. 
We  hope  that  the  present  work  might  help  in  understanding 
the  connection  between  the  two  perspectives.  This  is  another 
subject  for  future  research. 
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APPENDIX  A 


using  ip(t)  rather  than  <1/ it).  Thus,  according  to  the  procedure 
adopted  in  this  paper,  we  use  the  analytical  form  of  Eq.  (26). 
For  the  reader’s  convenience,  we  rewrite  this  expression 
here, 

=  (Al) 


It  is  important  to  stress  that  in  the  text  we  made  the  choice 
of  Eq.  (20)  for  <//(t),  in  this  way  determining,  through  Eq. 
(24),  a  form  for  that  departs  from  Eq.  (Al).  However, 
both  choices  lead  to  the  same  time  asymptotic  behavior  of 
Eq.  (14)  for  both  the  waiting  time  distributions. 

We  note  that  for  we  get  a  finite  value  <t>(°°)  =  (/r 

-1  )/T  corresponding  to  ip{ 0)  in  t  space.  As  in  Sec.  II,  we 
separate  the  kernel  into  two  contributions:  <L(f)  =  [(/r 
-  l)/T]<5(f)-i-<l>fl(f),  and  insert  the  Laplace  transform  of  sepa¬ 
ration  into  Eq.  (8),  to  write 

-  uip{ii)  -  {/i-  1)  +  (n~  \)ip{u) 

<vaw  = - i - 

1  -  (/<«) 

_  fe(«)  +  (^-  1  )*/<«)  ,A2, 

1  -  Lp(u) 


We  introduce  ipD(u)  as  the  Laplace  transform  of  the  distribu¬ 
tion’s  derivative  ip' (t)  =  -/i(/jL-\)Ttl~l / (r+f)^+1;  using  the 
Laplace  transform  of  an  inverse  power  law  and  substituting  it 
into  Eq.  (A2)  we  have 

_  /x)u(eu  -  E“) 

U1-* -  {/jl -\)T{\-/i){eu -E^)’ 

where,  as  usual,  for  simplicity  we  have  set  T=  1 .  Anticipating 
the  convolution  form  of  the  solution  we  cross-multiply  to 
obtain 


W  ^  —  (M  —  1)T(1  —  /*)(<?“  -  £"_,)]<!>„(«)  =  {/i-  l)[(/r, 

-  1)T(1  -  /i){eu  -  £^_i)]  -  /rT(-  /i)u{eu  -  £") .  (A4) 
Using  the  relation  [2] 

r°°  ta  TTQa 

I  - exp[-  ut\dt  =  - - (E‘a  -  eua),  a  >  -  1 , 

J  o  t  +  a  sin  ira 

(A5) 


and  setting  a=/i- 1,  we  construct 

R(t)=  r,  1  f  r-2$fl(f-f')rff'-— U-Dr(i 
r(/r  —  i)  Jo  77 

-/i)\  ~ — -<&„(f-  t')dt' ,  (A6) 

do  '  +  1 


In  this  appendix  we  show  how  to  obtain  the  order  of  the 
fractional  operator  in  the  GME  for  1  <  /i<  3 ,  using  the 
Laplace  transform  form  for  <I>(f)  given  by  Eq.  (8).  We  are 


where  R(t)  is  the  inverse  Laplace  transform  of  the  right  side 
of  Eq.  (A4).  Using  the  well  known  recurrence  relation  of  the 
T  function,  we  can  combine  terms  in  Eq.  (A6)  to  obtain 
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m  = 


M  1  - 


<&n(f-  t')dt' 


L-f„ 

1-1)  Jo 
i  r 

itItti** (A7) 


Tip. 

Let  us  consider  first  the  case  /jl>  2,  where,  in  the  limit  of  t 
S>1  (equivalent  to  t>T),  we  can  write  Eq.  (Al)  to  a  good 
approximation  as 


m  = 


-  f  t'^a(t-t')dt'  =R(t).  (A8) 

B(/r  —  1)  Jo 


Going  back  to  the  Laplace  transform  representation,  we  ob¬ 
tain  the  simplier  expression 

r(/t  — 2)<$a(n)  ,  ,/  „  u  a 

r(/«-i)  =<!*-»  r<i-K)(«  -V.) 

-Ma1-  1)F(- At)w(e“-£'^),  (A9) 
which  after  a  little  algebra  yields  for  (f>„(r/) 


^(m)  =  (//,-  l)(/r  —  2) 


(At- i)r(i  -/u,)— 


(A  10) 


Thus  using  the  inverse  Laplace  transforms  we  obtain  the  cor¬ 
responding  expression  in  the  time  representation: 


1 


1 


(1  +tY  (1  +  r)^1 
In  the  time  asymptotic  limit  we  get 


■  (All) 


(A  12) 


corresponding  to  a  fractional  operator  of  index  /3=/r- 2.  For 
the  sake  of  completeness,  we  also  give  the  expression  for  the 
GME  kernel  in  the  case  /i<2.  Proceeding  as  was  done  ear¬ 
lier  we  obtain 

*(»)  —  Si"  ^  +  (A  13) 

Finally,  we  want  to  point  out  that  the  expressions  we  are 
proposing  refer  to  values  of  /jl  which  are  not  integer.  We  are 
exploring  the  interval  [1,3].  Thus  the  the  expressions  we  are 
proposing  become  questionable  for  /j  =  2.  To  get  the  proper 
expression  for  /jl=2  we  have  to  study  expressions  like  those 
of  Eq.  (A2)  and  (A9)  at  fi=2  +  e,  do  a  Taylor  series  expan¬ 
sion  around  /jl= 2,  and  assign  to  ji-2  the  limiting  values 
reached  for  e— >0. 


APPENDIX  B 

This  appendix  is  devoted  to  establishing  the  accuracy  of 
Eq.  (53),  and  consequently  the  validity  of  the  assumption 


FIG.  4.  The  waiting  time  distribution  i/ft  ( t )  as  a  function  ot  time 
t.  The  dots  denote  the  exact  values,  and  the  lines  the  prediction  of 
Eq.  (55)  for  fi=2.3,T=  1.  Moving  from  bottom  to  top  in  the  right- 
hand  portion  of  the  figure,  fa=0.01,  1,  10,  100. 

that  the  beginning  of  the  first  laminar  region  overlapping 
with  t> 0  is  uniformly  distributed  between  t=—ta  and  t= 0. 
The  exact  expression  for  i//(  (f)  is 

f'a 

t/r,  (t)  =  I  dxG(ta  —  x)ifi(t  +  x) ,  (Bl) 

"  4o 


where 


G(t)  =  8(t)  +  t/*(f)  + 


-  n) ' 


S' 


2  [  dTf  i/A  ti  )  ■  f 

11=2  JO  J  Tj 


drn_\ilAt~  t„_i). 


dT2ifAr2 


(B2) 


It  is  straightforward  to  find  the  Laplace  transform  of  G(t). 
This  is  given  by 


G(«)  =  E  fa)”  = - (B3) 

n=o  1  -  ijAu) 

Thus,  the  Laplace  transform  of  Eq.  (Bl)  with  respect  to  ta 
reads 


(t)  : 


1 


l  -  t(Asa) 


<Ksa)  ■ 


f  e  S“yi/Ay)dy 

do 


.  (B4) 


By  numerically  anti-Laplace  transforming  Eq.  (B4)  we 
evaluate  the  time  dependence  of  the  exact  waiting  time  dis¬ 
tribution  of  age  ta,  Eq.  (Bl). 

In  Fig.  4  we  compare  the  exact  prediction,  evaluated  nu¬ 
merically,  to  the  heuristic  expression  of  Eq.  (53).  We  find 
that  at  small  and  large  values  of  ta  these  two  expressions 
coincide.  In  the  intermediate  region  they  do  not.  Neverthe¬ 
less,  we  think  the  agreement  between  the  two  expressions  is 
satisfactory  enough  for  the  purpose  of  this  paper. 
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